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Abstract: We study relevant deformations of an M = 1 superconformal theory 
which is an exactly marginal deformation of U(N) M = 4 SUSY Yang-Mills. The 
resulting theory has a classical Higgs branch that is a complex deformation of the 
orbifold C 3 /Z n x Z n that is a non-compact Calabi-Yau space with isolated conifold 
singularities. At these singular points in moduli space the theory exhibits a duality 
cascade and flows to a confining theory with a mass gap. By exactly solving the corre- 
sponding holomorphic matrix model we compute the exact quantum superpotential 
generated at the end of the duality cascade and calculate precisely how quantum 
effects deform the classical moduli space by replacing the conifold singularities with 
three-cycles of finite size. Locally the structure is that of the deformed conifold, but 
the global geometry is different. This desingularized quantum deformed geometry is 
the moduli space of probe D3-branes at the end of a duality cascade realized on the 
worldvolume of (fractional) D3-branes placed at the isolated conifold singularities in 
the deformation of the orbifold C 3 /Z n x Z n with discrete torsion. 



1. Introduction 



Gauge theories with Af = 1 supersymmetry (SUSY) in four dimensions are known 
to exhibit rich and beautiful dynamics. The phenomenon of Seiberg duality plays 
a central role in much of this [1]. One of the most striking examples of the rich 
physics of Af = 1 SUSY theories is provided by the duality cascade phenomenon 
uncovered by Klebanov and Strassler [2] . They investigated an Af = 1 gauge theory 
with SU (M + N) x SU (N) gauge group and matter fields in the bi-fundamental 
representation which, upon renormalization group (RG) flow, undergoes a succession 
of Seiberg dualities or a duality cascade. The number of colors in each factor drops 
repeatedly by M units until finally in the infrared (IR) we are left with a pure 
Af = 1 SUSY Yang-Mills theory that confines and produces a mass gap. 1 In the 
ultraviolet (UV), however, these theories are not asymptotically free and hence are 
not well-defined quantum field theories at all energy scales. 

In this article, we study the phenomenon of a duality cascade within a four 
dimensional UV-finite Af = 1 SUSY field theory. Our study was motivated by the 
program initiated by Berenstein in [4-6]. The theory in question is obtained from 
Af = 4 SUSY Yang-Mills theory with U (N) gauge group by deformations involving 
exactly marginal and relevant operators. Thus a UV-completion of the Klebanov- 
Strassler cascade can be achieved within a well-defined quantum field theory with 
gauge group of finite rank. Indeed, the theory we study has discrete branches where it 
flows, after a duality cascade, to an IR Af = 1 gauge theory with a mass gap. This IR 
theory (more precisely, its holomorphic sector) is identical to the theory studied in [7] 
which is quite similar to Af = 1* theory [8,9]. In the D-brane language, the theories 
that we study are realized on the worldvolume of D3-branes placed at isolated conifold 
singularities of a Calabi-Yau threefold which, in this case, is a complex deformation of 
the orbifold C 3 /Z n x Z n . The endpoint of the duality cascade in the field theory has 
a geometric dual description where the D-branes have disappeared and the Calabi- 
Yau geometry is deformed and desingularized by the appearance of three-cycles of 
finite size. One of our aims is to obtain the precise algebraic description of this 
quantum deformed geometry with deformation parameters determined exactly in 
terms of the parameters of the gauge theory. Although near the singularities its 
local structure is that of the deformed conifold, the global structure is very different. 
From this deformed geometry we will also understand how the theory of [2] emerges 
as a limiting case when one of the finite three-cycles in the geometry is sent to 
infinity and we are left with a deformed conifold geometry. The main tool that we 
will employ to understand this deformed geometry (which is also the moduli space 
of probe D-branes at the end of the cascade) and the gauge theory superpotential, 

1 See the recent work [3] for further discussions on the IR universality class of these theories. 
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is the holomorphic matrix model approach of Dijkgraaf and Vafa [10-12]. 



The theories of interest to us are obtained from the so-called u (3- deformation" of 
the M = 4 theory with U(N) gauge group. The /9-deformation constitutes one of the 
two J\f = 1 supersymmetric, exactly marginal deformations of the M = 4 RG fixed 
point identified by Leigh and Strassler in [13]. It leads to an M = 1 superconformal 
field theory (SCFT) with three adjoint chiral supermultiplets X, Y and Z which 
obtain the tree level superpotential interaction 2 

W = XTr[XYZ-qXZY]; q:=e ifi . (1.1) 

where A and q are complex bare parameters. We also introduce the complex bare 
gauge coupling of this theory r = Ani/ ' g\ M + 9/2tt. We will always restrict attention 
to the case where q is a primitive n-th root of unity, q = exp(2m/n). In order for 
the theory to be conformally invariant, the other coupling A has to be some fixed 
(but unknown) function of (3 and the gauge coupling [13] such that as we approach 
the M = 4 fixed line (f3 — >• 0), we also have A — > 1. The gauge-coupling r is a free 
parameter of the theory and it was shown in [7] that the theory exhibits electric- 
magnetic duality which acts as SL(2, Z) transformations on the gauge coupling: 3 

ct + a ct + a 



Dynamical phenomena of interest ensue when we break conformal invariance 
by perturbing the above Af = 1 SCFT, with q = exp(2ni/n), by certain relevant 
operators. Specifically, we will consider adding to the tree level superpotential, a 
mass term for one of the fields of the following form 

W = XTi[XYZ -e^XZY -e(\X 2 -aX)}; n>2. (1.3) 

More complicated polynomial deformations can also be considered leading to similar 
IR dynamics, but these spoil the UV-behaviour. We also consider slightly more 
general renormalizable deformation in the penultimate section. 

Let us now explain the essential ideas underlying the dynamics of the above 
theories. It is well known that when q is an n-th root of unity the /9-deformation 
of the U(N), M = 4 theory can be realized in IIB string theory as the worldvolume 
theory of iV D3-branes placed at a C 3 /Z n x Z n orbifold singularity with a single 
unit of discrete torsion [14,15]. The moduli space of vacua of the /9-deformed theory 

2 We work with a normalization where the kinetic terms of the chiral multiplets do not have a 
factor of 1 / 9ym m front- 

3 Actually, in order to write the transformations in this way the gauge coupling undergoes an 
algebraic renormalization: we refer to [7] for details. 
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contains Coulomb and Higgs (and mixed) branches. It follows from the F-term 
equations that, depending on the rank N of the gauge group, the Higgs branches 
within the moduli space are symmetric products of a certain number of copies of 
C 3 /Z n x Z n . (For example, when N = mn there exists a Higgs branch where the 
gauge group is broken to a U(m) subgroup and the Higgs branch is a symmetric 
product of m copies of C 3 /Z n x Z n ). The orbifold has three singular (complex) 
lines which are fixed under the action of the orbifold group, and along which the 
Coulomb and Higgs branches of the D-brane gauge theory meet. The separations of 
the D-branes along the fixed lines are the Coulomb branch moduli. At these singular 
fixed lines there are fractional D-branes bound to the singularities 4 , which can then 
annihilate in groups of n to form a whole D-brane free to move off into the orbifold 
bulk (the Higgs branch.) 

Once the theory has been deformed by relevant operators, the fixed line singu- 
larities above can be smoothed out by complex structure deformations. From the 
point of view of the gauge theory on the D-branes this corresponds to deforming 
the superpotential. Deformations that resolve the fixed lines leaving behind only 
isolated conifold singularities were first discussed in [15] . The deformation ( |1.3|) that 
we consider was discussed in [6] and [16,17]. After a suitable field redefinition it is 
equivalent to introducing masses for all three multiplets, a deformation that has been 
studied for general complex (3 in [7]. Importantly, for the theories with (3 = 27r/n, 
this class of deformations does not lift the Higgs branch moduli space, rather it leads 
to a moduli space which is a complex deformation of the C 3 /Z n x Z n orbifold. Clas- 
sically, the effect of the deformation is to smooth out two of the three fixed lines 
into a single (complex) hyperboloid, while the third fixed line is completely resolved 
leaving behind [n/2] isolated conifold singularities. 5 ^From the string theory point of 
view, once we have isolated conifold singularities in the spacetime geometry we can 
imagine having fractional branes stuck at these points. Indeed, in the field theory 
perspective, a bulk D-brane is characterized by an n-dimensional irreducible repre- 
sentation of the algebra following from the F-term conditions implied by ( |1.3j ). These 
representations become reducible precisely at each of the [n/2] points in moduli space 
corresponding to the conifold singularities. At any one of the conifold singularities 
labelled by an integer p, a bulk D-brane splits up into two lower dimensional irre- 
ducible representations of dimension p and n — p respectively. These are the two 
types of fractional D-branes bound to the singularity. 

The appearance of isolated conifold singularities and fractional branes pinned 
at them naturally points towards the possibility of realizing M = 1 gauge theory 
dynamics of the kind studied by Klebanov and Strassler [2] . This was already realized 

4 See [15] for a discussion on the interpretation of these fractional branes. 
5 Note that [x] is the largest integer less than or equal to x. 
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by Berenstein in [4] and [6]. In this paper we will be able to make this precise. First 
we will show that for the U(N) gauge theory with classical superpotential ( |1.3|) , the 
low energy effective theory at the isolated singular points is a U(N + ) x U(N~) gauge 
theory (here N = pN + + (n — p)N~) with bi- fundamental matter and a quartic 
superpotential. In the D-brane language N^ 1 is the number of fractional branes of 
each type. In an appropriate regime of parameter space which results in a separation 
of scales, this theory exhibits a duality cascade that terminates in an Af = 1 gauge 
theory with a mass gap. It is important to note that the low-energy theory is not 
pure M = 1 SUSY Yang-Mills theory. In fact it has additional massive fields and 
is an M = 1* type theory which was studied in [7]. In particular the low energy 
theory develops a superpotential and gaugino condensate both of which are non- 
trivial modular functions of the microscopic gauge coupling r, the same functions 
that were encountered in [7]. Only in a specific decoupling limit do we obtain pure 
M = 1 SUSY Yang-Mills. All this is encoded rather beautifully in the geometry of 
the quantum deformed moduli space which we are able to obtain exactly. 

The classical moduli space has conifold singularities where fractional branes are 
pinned. From general considerations [2,18,19] in the quantum theory we expect the 
fractional branes to deform the geometry further. At the endpoint of the cascade 
we expect that the (fractional) D-branes have disappeared, the conifold singularities 
have been replaced by >S 3 's of finite size resulting in a quantum deformed moduli space 
seen by a probe D-branes at the endpoint of the cascade. Since all the holomorphic 
data of an Af = 1 gauge theory is encoded in a corresponding matrix model [10-12], 
we should be able to determine this quantum deformed moduli space using the matrix 
model approach. We achieve this using the fact that the planar sector of the matrix 
model for our theory is in fact exactly solvable [7,20-22]. Remarkably, in the matrix 
model one can understand the duality cascade of the field theory in rather simple 
terms. It turns out that the matrix model description of the field theory vacuum of 
interest with U(N + ) x U(N~) gauge symmetry, only depends on a single parameter 
corresponding to the glueball superfield of the U(\N + — N~\) gauge group factor 
appearing at the end of the cascade. Here \N + — N~\ has the interpretation as the 
total number of fractional branes pinned at the conifold singularity. This particular 
feature was also encountered in [23] in the context of the matrix model description 
of A\ quiver theory of [2]. 

For us the deformed Calabi-Yau geometry emerges as a fibration over a Riemann 
surface that is completely determined by the matrix model. This quantum geometry 
is a deformation of the original classical Calabi-Yau space with conifold singularities, 
and we find that the quantum deformation parameters can be determined exactly as 
a function of the microscopic gauge coupling r. The conifold singularities are now 
replaced by three-cycles whose sizes are exactly calculable, nontrivial functions of 
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the gauge coupling or equivalently the glueball superfield. An important feature of 
these deformed geometries is that each of the "blown-up" S 3, s also has a finite dual 
three-cycle. It is precisely in the limit that this dual cycle is sent to infinity that we 
are left with just the usual deformed conifold geometry which is dual to pure Af = 1 
SUSY Yang-Mills theory 

The layout of this paper is as follows. In Section 2, we review the classical 
vacuum structure of the field theory and its classical moduli space which has isolated 
singularities. We describe the classical solutions that correspond to different types of 
fractional branes stuck at singularities. In Section 3, we obtain the classical effective 
action for the light modes at the singular points in moduli space and observe that it 
corresponds (at weak coupling) to the starting point of a duality cascade. In Section 
4, by solving the Dijkgraaf-Vafa matrix model associated to this theory we find the 
quantum deformed moduli space where the deformation parameters are determined 
precisely in terms of parameters of the gauge theory. Section 5, considers theories 
with a slightly more general renormalizable deformation including the terms TrXY 
and Tr XZ. These theories can easily be solved along the lines we set out in earlier 
sections. Finally, in Section 6 we conclude with some remarks on directions requiring 
further study. In Appendix A, we derive the quadratic loop equations for the matrix 
models in question and point out a way to obtain the general loop equations. Various 
useful results involving elliptic functions are listed in Appendix B. 

2. Classical Vacuum Structure 

In this section we address the question of the classical vacuum structure of the (3- 
deformed theory with superpotential (|1.3|). In fact, following [6] we will consider a 
generalized version of the deformation in (|1.3| ). Since we are primarily interested in 
the holomorphic structure we shall follow the rule of solving the F-flatness conditions 
modulo complex gauge transformations. Subject to certain well-known caveats, such 
solutions are then equivalent to the solution of both the F- and D-flatness conditions 
up to a complex gauge transformation. We consider the gauge theory perturbed by 
a polynomial superpotential 6 

I 

W = Tr[XYZ - qXZY + V(X)\ , V{x) = ~ xi ■ ( 2 - X ) 

jj^O mod n ^ 

The fact that we exclude terms in the superpotential proportional to X\ j = mod n 
is the restriction of [6] and ensures that the deformation does not lift the Higgs 

6 In the rest of the paper we set A = 1. It could easily be reinstated by by a re-scaling of the 
fields. 
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branch. Note that only superpotentials for which I = 2 are actually renormalizable 
and we shall mostly have this case in mind because we want theories which flow to 
the superconformal fixed-point in the UV. Note also that in this case once we remove 
the linear term in X by a shift, the superpotential is equivalent to 



W = Tr 



XYZ - qXZY - Ql(1 q) YZ + ^X 2 } , (2.2) 

0,2 2 



which with a trivial rescaling is precisely the theory investigated in [7]. However, for 
illustrative purposes it is also useful to consider the q = — 1 theory with £ = 3 where 
the superpotential is actually marginally irrelevant, or non-renormalizable. 

These perturbed theories and their moduli spaces of vacua have been studied 
in [5,6,16,17]. In what follows, we shall draw heavily on this analysis which shows that 
the classical moduli space of vacua on a Higgs branch of these theories corresponds 
to be a 3-complex dimensional Calabi-Yau space X which is a complex deformation 
of the orbifold with discrete torsion (i.e. a deformation of C 3 /Z n x Z n ). We will 
further see that the low-energy description of the gauge theory involves both discrete 
vacua and a moduli space with several branches. On some of these branches, the 
gauge theory exhibits a subsequent renormalization group cascade of the Klebanov- 
Strassler type [2], and on others the theory flows to the superconformal fixed point 
of Klebanov and Witten [24] . In addition there are isolated vacua where the theory 
confines without undergoing a duality cascade. It is worth mentioning that the 
existence of this rich vacuum structure is then inherited in the theories related by 
the SL(2, Z) group. In particular under S-duality, the Higgs branches are mapped 
to confining branches of the dual theory with modified gauge coupling 1/r and a 
deformation parameter 27r/(rar) [7]. 



2.1 The case q = — 1 with cubic deformation 



It is instructive to first specialize to the simplest example, namely the theory with 
n = 2 (and q — — 1) and subsequently describe the situation for q a generic n-th 
root of unity. Consider then the superconformal theory with q = — 1 in Eq. ( |1.1| ) 
perturbed by a tree level superpotential for one of the fields X. We restrict attention 
to a perturbation of the form 

V(x) = -e(x 3 /3-a 2 x). (2.3) 

This perturbation satisfies the restriction implied in Eq. ( p.l| ) which only requires 
that it should be a polynomial in odd powers of X. As we discussed earlier, a 
cubic superpotential would be marginally irrelevant once quantum effects are taken 
into account. This fact renders this example less interesting since it spoils the UV 
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behaviour. However, it is useful to discuss it because of its simplicity and because 
it exhibits all the interesting infrared dynamics that we are after. Later we shall 
consider theories with n > 2 with quadratic superpotentials which exhibit the same 
properties but which are also renormalizable. We shall also consider a more general 
class of theories, including q — — 1, which are renormalizable. 

The F-term equations for the theory with q = — 1 are 

{X,Y} = {X,Z} = 0, {Y, Z} = e{X 2 - a 2 ) . (2.4) 

The classical moduli space of this theory can now be obtained following the methods 
discussed in [5,6, 16]. First of all, complex gauge transformations can be used to 
diagonalize, say, X. This leaves unfixed diagonal gauge transformations as well as 
those generated by the Weyl group which act by permutation of the two diagonal 
elements. We will have more to say about these transformations shortly. For the 
above choice of potential, the F-term equations are solved by the following set of 
fields 

X = xa 3 ; Y = ya 2 - Z = z a 1 + z' a 2 (2.5) 
where the a 1 are Pauli matrices with 

2 yz' = e (x 2 -a 2 ). (2.6) 

However, this solution is fixed by a Z 2 x Z 2 subgroup of the gauge group (the product 
of the Weyl group and a subgroup of the diagonal transformations). Explicitly the 
two generators are 

9i = a 3 , 6 2 = a l . (2.7) 
Modding out by these transformations is equivalent to identifying 

01 : (x, y, z) ~ (x, -y, -z) , 9 2 : (x, y, z) ~ (-x, -y, z) . (2.8) 

In the D-brane interpretation of [16, 17] this irreducible 2x2 dimensional rep- 
resentation should be thought of as parameterizing the moduli space of a single 
D3-brane probe moving on a 3-complex dimensional Calabi-Yau space X. In order 
to find an explicit form for this classical moduli space we need to introduce gauge 
invariant variables which automatically incorporate the identifications (|2.8|) . This 
can be done by first enumerating the set of elements that commute with all the gen- 
erators of the algebra in Eq. ( |2.4| ). One finds X 2 , Y 2 and Z 2 . These 3 generators 
are therefore proportional to the identity matrix and we can introduce the gauge 
invariant variables (u, v, w) via 

u ■ l[2]xp] — X 2 , v ■ l[2]x[2] — Y 2 , w ■ 1[2] X [2] = Z 2 . (2-9) 
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Clearly 




w = z 2 + z' 



a 



(2.10) 



There is an additional central element: 



t-l 



[2]x[2] — 



XYZ + \X V'(X) . 



(2.11) 



The four variables (t,u,v,w) are not independent, rather they satisfy the relation 



This algebraic equation describes a (non-compact) 3-complex dimensional Calabi- 
Yau space X in C 4 which is precisely a complex deformation of the C 3 /Z2 x Z 2 , 
the orbifold with discrete torsion. Indeed, when the deformation vanishes we have 
t 2 = —uvw which is the well-known algebraic description of the orbifold. The un- 
deformed orbifold with discrete torsion has (classical) singularities corresponding to 
the (complex) fixed lines where any two of u, v and w vanish. These fixed lines 
correspond to solutions (|2.5|) (with e = 0) fixed by 9\, 9 2 or 9\9 2 , so y = z = 0, 
x = y = and x = z = 0, respectively. 

The addition of the superpotential perturbation (e 7^ 0) resolves the orbifold 
singularities, but only partially. In particular, it is easy to see that the singular 
lines associated to the undeformed solution fixed by the gauge transformations 9 2 
and 9\9 2 (u = v = and u = w = 0, respectively) merge. For instance, the 
singular line associated to 9 X 9 2 becomes x = z = and yz' = —ea 2 /2, or u = and 
wv = (— ea 2 /2) 2 . However, the third singular line associated to the solution fixed 
by 9\ is reduced to an isolated singular point located at {u = a 2 ,v = 0,w = 0). 
This singularity is locally of the conifold type. 7 In fact precisely as we approach this 
singular point from the bulk, the two-dimensional solutions to the F-term equations 
above become reducible in the sense that there are now two one-dimensional solutions 
y = z = and x = ±a. In the brane interpretation, the fact that the two-dimensional 
irreducible representations of the deformed algebra split up into two one- dimensional 
representations can be interpreted as the fractionation of a bulk brane at the conifold 
singularity. In particular each one-dimensional solution corresponds to a fractional 
brane stuck at the conifold singularity while the irreducible two-dimensional solution 
describes the possibility for the two types of fractional branes at the singularity to 
combine into a whole brane and move away from the singular point into the bulk of 
the deformed geometry which is also the Higgs branch of the D-brane gauge theory. 
(Note that the one- dimensional solutions only exist at u = a 2 and so an individual 
fractional brane can only be pinned at this singularity.) In the superconformal theory 

7 This can be seen by introducing the new variable 5 = u — a 2 and focusing attention on S, v, w <ti 



uvw + —ru{u — a 2 ) 2 . 



(2.12) 
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with e = this conifold singularity becomes a point on the fixed line v — w — 
where the Coulomb branch meets the Higgs branch. 



2.2 The general case n > 2 

The picture that emerges for generic n > 2 is quite similar. If we are only interested 
in the low-energy dynamics we can always allow arbitrary polynomial perturbations 
satisfying the requirement (|2.1|) resulting in classically deformed moduli spaces that 
will be described below. The restriction to renormalizable terms will become impor- 
tant in the context of what we have to say in the next section. 

The F-term equations for the perturbed theory are 



To start with there are two Coulomb branches where X = Z = and Y diagonal and 
X = Y = and Z diagonal, respectively. We will not be interested in these branches. 
The remaining solutions can be built out of irreducible blocks of size up to n. Blocks 
of dimension p < n come in sets {9^ } while the block of dimension n is unique. 
In the q = — 1 case above these were the two one-dimensional solutions and the two- 
dimensional irreducible solution, respectively. The solutions with dimension p < n, 
the {9^}, will describe fractional branes pinned at conifold singularities. The other 
irreducible n-dimensional solution D^™-* will describe whole branes that are free to 
move off the singularities and explore the bulk of the classical moduli space. 

In the following a central role will be played by the function f(x) defined in 
terms of the potential by 



XY - qYX = , ZX- qXZ = , YZ - qZY = -V\X) 



(2.13) 



xV\x) 



/(*) 



f(xq) . 



(2.14) 



Note that 




jj^O mod 



n 



(2.15) 



and has the same order as V(x). The two classes of solutions to ( p,13|) are more 



specifically 
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i) 9\\ i « = d^ p \ of dimension p < n: 



n 



where ^f 1 is a root of the polynomial equation 



\ 



(p) 



1 
1 



q p-lJ 



/0 *i 

z 2 



(p) 



1 oj 



f(x) = f(x<f) 



z p -i 

o / 

(2.16) 



(2.17) 



excluding x — 0. 



(ii) = x n ,f = y n ,w = z n ], of dimension n: 



/I 



\ 



V 



,n-l 



y(n) 



1 
1 

V 



1\ 



z (n) 



1 0/ 



/o Ci 

o C 2 



Cn-1 

o / 

(2.18) 



with rij=i — 1- The actual values of the Q are fixed but will not be required. 

Since x, y and z are not fixed the representation 9t( n )[u, v, w] has a moduli space 
which can be obtained, as in the q — — 1 case, by determining the relation between 
the elements of the center of the algebra, namely 



U ' l[n]x[n] — X n , V ■ l[„]x[„] — Y n , W ■ l[ n ]x[n] ~ 5 ^ ' l[n]x[n] 

The set of variables (t, u, v, w) satisfy the relation 

n 

uvw = (-lf +1 ]J{t- f(xq*)) = (-l) n+1 F(t, u) 

3=1 



XYZ + /(X) . 

(2.19) 



(2.20) 



Notice that the right hand side of this equation is actually invariant under x — > qx 
and can therefore be only a polynomial F(t,u) in the invariant u = x n . This is 
a 3-complex dimensional Calabi-Yau space X which is a deformation of the orb- 
ifold C 3 /Z n x Z n with discrete torsion described by the algebraic equation t n = 
(-l) n+1 uvw. 

As in the q — — 1 case the singular fixed lines of the orbifold are partially re- 
solved. Rather than investigating these singularities by considering the action of 



10 



gauge transformation as we did for q — — 1, the structure of the singularities can 
equivalently be determined directly from the Calabi-Yau geometry (|2.20|) , as in [6], 
by demanding that the partial derivatives with respect to u, v, w and t vanish. The 
first kind correspond to u = along with 

vw = d u F(t,u)\ u=0 , (2.21) 

which is a deformation of the two fixed lines u = v = and u = w = 0. The 
other class of singularities are the ones that interest us. They are isolated points in 
the geometry where v = w = and the the n-sheeted algebraic curve, or Riemann 
surface, £ defined by 

F(t,u) = (2.22) 
becomes singular. These points are simply given by the roots of 

f(x) = f(xq p ) (2.23) 

for any p < n, excluding x = 0. Hence the isolated singular points are precisely at 
the locations of the fractional branes associated to the representation 9\[ p \ It 

Earlier we denoted these roots e\ , where i counts the degeneracy: i — 1, . . . , S p \ 
The fractional branes are pinned at a conifold-type singularity at 



u = u- 



(p) _ (Ap) 



(ef) n , v = w = 0. (2.24) 



Let us enumerate the number of these isolated singular points. Recall that we are 
counting the roots of / '(x) — / (xq p ) which is a polynomial equation of degree £, except 
when £p = mod n, in which case it has degree I — 1. So the number of roots, apart 
from x = 0, is 

I £ — 2 £p = mod n 
d (p) = { „ (2-25) 
- 1 otherwise . 



Note also that there is a symmetry 



t~ P) = (2-26) 



and d( pS> = S n ~ p \ (These equivalences are valid when p = n/2 for n even.) This 
means that the fractio 
point in the geometry: 



means that the fractional branes associated to 91-^ and , n . lie at the same 



^ = ■ (2-27) 

This is very familiar from the brane interpretation of the basic conifold theory: frac- 
tional branes come in pairs. In order to avoid this cumbersome notation, we shall 
label the singularities, or pairs of mirror fractional branes, by the index /x = (p,i). 
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In particular, for a given \x the associated value of p is 1 < p^ < [n/2] and we denote 

{n ~ p) as<R+ 
dW+i-i as m 



the mirror pair 9\ {p) and IK^,. as and 9tT. Finally 



^ = ^ = u ^ 1 _., ^ = 1,...^. (2.28) 
Note that the number of singularities, or pairs of fractional branes, is equal to 

, (n/2) Kn-l)/2] 
3=1 



However, can we see that if we have two fractional branes, one of each type, 
then it can become a "complete" brane and move off into the bulk? The answer 
is yes because at (u = u^,v — 0, w — 0) the representation ^ n '[u,v,w) is actually 
reducible: 

9tf»)[ M = u^v = 0,w = 0] ~ ^+ ©1K M . (2.30) 

Let us consider in more detail the relation of these singularities to the algebraic 
curve S defined by F(t,u) = 0. We can think of £ as an ra-fold cover the u plane. 
There is a n th -order branch cut running from u = to u = oo which joins all the 
sheets. In addition for u = a pair of sheets touch at a point. So as one expects 
the curve £ is singular at these points. In the quantum theory the picture emerging 
will be that the presence of a net number fractional branes, i.e. the number of 
representations fR^ minus the number of has the effect of deforming the isolated 
conifold singularity in X which at the level of the curve £ involves the opening up 
of the point u = u M into a cut joining the pair of sheets. 

Now that we have described the irreducible blocks, we can now write down a 
general classical vacuum. If one starts with a U(N) gauge group then a general 
vacuum is associated to a reducible representation 

©£ =1 (iV;-^©Ar -9t;) , (2.31) 

where 

N = E (p^ + ( n - p/W) • ( 2 - 32 ) 

M=l 

Notice that we didn't include the block explicitly because the vacuum ( |2.31| ) 
lies at a point in a large moduli space where up to |iV+ — N~\, for each pair 9^ of 
fractional branes, can move of into the bulk and explore the space X. 
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3. The Classical Effective Action and RG Cascade 



We have seen that certain deformations of the /9-deformed M = 4 theory lead to 
deformed (classical) moduli spaces for D-branes with isolated conifold singularities 
where fractional branes can be stuck. We will now ask what the effective gauge 
theory is on stacks of branes some of which remain stuck at the isolated points while 
others have moved off to explore the full moduli space. As in the last Section it is 
useful to consider the q = — 1 case first. 



3.1 The case q — — 1 with cubic deformation 

Consider the U(N) gauge theory with the cubic plus linear superpotential 
Recall that there is a single isolated singularity in X at u = a 2 , v = w = 0. The two 
fractional branes are distinguished by x = ±a. So let us consider a vacuum with N + 
(fractional) branes with x = +a and N~ branes with x = —a (fractional branes of 
the other type) with N + + N~ = N. In this vacuum 

x= /^-Wi^] o \ Y = z = (31) 

V a ' l[iV~]x[JV~]/ 

Note that as we remove the deformation smoothly V — > 0, these points lie at the 
intersection of the Coulomb and Higgs branches of the parent superconformal theory. 



The classical vacuum (|3.1|) preserves a U(N + ) x U(N ) gauge symmetry. 



The low-energy classical U(N + ) x U (N ) gauge theory at these points in moduli 
space is obtained by expanding in fluctuations about the above solution, 

IX = ft ° V SY= ft • «? = (l 1 i 2 ) • (3-2) 



where X±, Yi and Z\ transform in the adjoint representation of the unbroken U(N + ), 
while X 2 , Y 2 and Z 2 are in the adjoint of the U(N~). In addition the fields A\ and 
Bi fill out four bi-fundamental multiplets transforming as (N + ,N~) and (N + ,N~) 
respectively. 

The nature of the low-energy theory can now be easily determined by expanding 
out the classical superpotential in terms of the fluctuations and we find 

W = Tr[X 1 (A 1 B 2 + A 2 B X ) - taX\ - \x\ + 2aY 1 Z 1 + X 1 (Z 1 Y 1 + Y X Z X )\ 

e ( 3 - 3 ) 
+ Tr[X 2 (B 1 A 2 + B 2 Ai) + taX\ - -X\ - 2aY 2 Z 2 + X 2 (Y 2 Z 2 + Z 2 Y 2 )} . 
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Note that the fluctuations Y i2 and Z lj2 have a mass set by the scale a while those 
associated to X\ and X 2 have a mass set by the scale ea. Assuming e C 1 and the 
gauge coupling is weak enough, it makes sense we can integrate out Y\ )2 and Z^ 2 at 
the classical level to arrive at the effective superpotential 

W eS = Ti[X 1 (A 1 B 2 + A 2 B X ) - eaXl\ + Tr^B^ + B 2 A X ) + eaX\\ , (3.4) 

where we have neglected the cubic coupling of X 12 which will be irrelevant in the IR. 
The effective theory for the light modes is a mass deformation of the U(N + ) x U (N~) 
Af = 2 theory with equal and opposite sign masses for the two adjoint scalars. This 
is precisely the starting point of the Klebanov-Strassler duality cascade. We can also 
now integrate out X lj2 to arrive at 

W cS = -^Tr(A5^ fe 5*KV 7 . (3.5) 

For a fixed given value of iV we can vary iV + (or N~) and explore different vacua 
of the theory. In particular when the unbroken rank iV is even we can choose the 
vacuum with N + = N~ = N/2. Here the resulting U(N/2) x U(N/2) theory with 
bi-fundamentals and the quartic superpotential above flows to the superconformal 
conifold theory of Klebanov and Witten [24]. As pointed out in [24], the quartic 
superpotential is actually an exactly marginal operator at the IR fixed point. Thus 
we have found a flow that interpolates between the C 3 /Z 2 x Z 2 theory and the 
Klebanov- Witten conifold theory. 

For any given N = iV + + N~ we can sweep through distinct classical vacua by 
varying N + . In particular, when N + = 0, we simply have a U(N) gauge theory with 
Af = 1 SUSY with massive adjoint scalars. The SU (N) subgroup of the U (N) gauge 
group confines and generates a mass gap and in many ways is similar to the confining 
vacuum of Af = 1* gauge theory [8,27]. In an appropriate decoupling limit for the 
adjoint scalars, the theory reduces to pure Af = 1 SUSY gauge theory. 

The remaining vacua labelled by generic values of N + contain rich and com- 
plicated dynamics. In particular, at least at weak coupling g\ M <C 1 (the gauge- 
coupling of the parent q = — 1 SCFT), the low-energy U(N + ) x U(N~) theory with 
bi-fundamentals and a quartic superpotential shares all the features of the theory 
studied by Klebanov and Strassler [2]. Under renormalization group flow, the indi- 
vidual gauge group factors will undergo successive Seiberg dualities that terminate 
at an Af = 1 theory with a mass gap and gaugino condensation, with or without 
probe branes propagating on a deformed and desingularized moduli space. We will 
return to a more detailed quantum description of these vacua after demonstrating 
how they arise in the theories with n > 2. 
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3.2 The case n > 2 with quadratic deformation 



At the outset we remark that there is a key difference between the cascade theories 
encountered here (and above) and that of Klebanov-Strassler. The theory consid- 
ered in [2] has the property that one of the gauge group factors has a Landau pole in 
the ultraviolet. (One could imagine achieving a UV completion by a reverse duality 
cascade, but in that case the ranks of the gauge group factors would grow without 
bound in the UV 8 .) The theories we consider have the property that they are ob- 
tained as deformations of a four dimensional M = 1 superconformal field theory with 
U(N) gauge group and therefore any potential duality cascade is embedded within a 
well-defined UV-theory. In this context, note that although the deforming potential 
for q — — 1 was classically marginal (i.e. cubic), it is most likely to be rendered 
irrelevant by quantum effects. However, as we will see below for the theories with 
n > 2 the deforming potentials need only be quadratic — hence renormalizable — to 
obtain an RG cascade in the infrared. 

We consider the superconformal theory perturbed by the superpotential 

V(x) = -e{x 2 /2-ax). (3.6) 

This superpotential yields exactly one fractional brane solution for each irreducible 
representation of dimension p < n (except for p = n/2 when n is even) stuck at 

u = u(p) = ( a r^)"' ^ = ^ = (3.7) 

in X. We now consider the theory with a U(N), N = pN + + (n — p)N~ gauge 
group in a vacuum with N + fractional 9^ p ) branes and N~ fractional 9^™~ p ) branes. 
We are going to argue that as in the q — — 1 case the classical effective theory is a 
U(N + ) x U(N~) theory with massless bi-fundamentals. First-of-all, the gauge group 
is obviously broken to U(N + ) x U(N~). This is clear from the from of the VEV 
assigned to the adjoint field X: 

x= (l [N+]x[N+] ®X^ 

V o i [JV - ]x[JV - ] ®x("-p) 

But the question is whether there are the right number of massless bi-fundamental 
fields. Indeed the answer is yes, as we can see by first considering the simpler case of 
the theory realized on one brane and one brane. As we have already seen, 

these two fractional branes make a complete brane that can move off onto the Higgs 
branch and thus there must be additional massless moduli. These additional moduli 

8 For a discussion of the UV behaviour of various theories exhibiting the phenomenon of duality 

cascades, see [25, 26] 
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correspond to the non-zero elements of Y and Z in the representation 9^ n )[w, v , w] 
that are zero in 

fftfp) ©<h(»-p). More specifically if j is a matrix with a 1 in position 
(ij) and zeros elsewhere, then the elements in question are E p+ i tP and E^ n in Y and 
and E n>1 in Z. 

Therefore, when we have N + coincident Ot^ branes and N~ !JK( n_p ) branes at 
the singular point there will be 4 x N + x N" such massless fields which will fill out 
a bi- fundamental multiplet in the (N + ,N~) © (N + ,N~~) of the low-energy product 
gauge group. To be more specific the low-energy degrees of freedom can be encoded 
in the following set of fluctuating modes 

«=C, 1 ^M+te- (3.9c) 



B 2 ® E n _ pA 

Here X\ transforms in the adjoint representation of the unbroken U(N + ) while X 2 
transforms in the adjoint of U (N~). We have also defined = diag(l, q, . . . , q p ~ 1 )- 
The modes (A i: Bi) are the bi-fundamentals transforming in the (N + ,N~) and 
(N + ,N~) representations, respectively. The remaining fluctuations (modulo gauge 
transformations) are denoted 5X, etc.. The classical superpotential for the fluctua- 
tions at the chosen vacuum is then 

i- = TrpC^^ - A 2 B 1 )) + Tr(X 2 (B 1 A 2 - B^)) 

S • 2 P (3.10) 

+ € -.t«(TrXl-TrXl) + ~5W. 
I 1 — q z 

The remaining factor 5W includes all the additional fluctuating modes. Importantly 
all these modes gain a mass at a scale a and so for sufficiently large a compared with 
e there is a separation-of-scales between these tilded modes and the ones we have 
singled out and we can ignore them — at least at sufficiently weak coupling. Note 
that the effective couplings of the U(N + ) and U(N~) factors are pr and (n —p)r, 
respectively. 

We might have expected the above picture to be true since from a geometric 
point-of-view taking a ^> e zeros-in on the singularity which is then well approx- 
imated by the deformed conifold and the usual story of the duality cascade and 
deformed conifold should then be applicable [2]. Indeed upon integrating out the 
massive adjoints we obtain precisely, the quartic superpotential 

W eS = -1 • l^-TriABjAkB^ 1 . (3.11) 
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Note that we have obtained what will be the starting point for a duality cascade under 
renormalization group flow, specifically when the gauge coupling of the ultraviolet 
theory is taken to be sufficiently weak in order that the effective superpotential ( ft.lOQ 
is valid over a certain range of energies between the decoupling of the additional 
modes and the onset of the cascade. What is new and remarkable is that the duality 
cascade has been embedded in a four-dimensional theory with a perfectly well-defined 
UV limit given by the superconformal fixed point. The RG flow away from the fixed 
point is defined by the mass deformation ( |3.6|) . The end-point of the cascade is a 
U(\N + — N~\) theory which confines. As noted earlier, when N + = N~ the theory 
flows to the conifold superconformal theory realizing a flow between the C 3 /Z n x Z n 
orbifold and the Klebanov-Witten SCFT. 

It is clearly possible to generalize the discussion to a vacuum with more than 
one type of fractional brane pair. In this case, at low energy the theories associated 
to each pair will decouple from each other. 



4. Engineering the Quantum Geometry from the Matrix Model 

We have seen how the classical low-energy description of the perturbed SCFT leads 
to a picture of the gauge theory being realized on (fractional) D-branes placed at 
conifold singularities in a geometry X which is a deformation of the C 3 /Z n x Z n orb- 
ifold. This deformed geometry was precisely the Higgs branch moduli space of the 
gauge theory. ^From the classical picture we also concluded that at generic classical 
vacua the theory undergoes strong-coupling quantum dynamics such as successive 
Seiberg dualities leading to a duality cascade and/or confinement with gaugino con- 
densation. We expect that the nonperturbative quantum effects (i.e. the presence 
of fractional branes) should result in a further deformation of the geometry above 
which desingularizes the isolated conifold points. This was exactly the picture at 
the endpoint of the RG cascade in [2] wherein the associated geometry was just the 
deformed conifold. But this cannot be the whole story in the class of theories that 
we are studying as they are actually UV finite theories that confine in the IR with 
or without an intermediate duality cascade. In a certain decoupling limit for the 
adjoint scalars they will reduce precisely to the theory of [2]. Correspondingly, the 
deformed Calabi-Yau geometry X (which will be the moduli space seen by probe 
branes at the end of a cascade) will only reduce to the deformed conifold in a limit 
where some 3-cycle is sent to infinity. 

^From the work of Dijkgraaf and Vafa [10-12], it is now well-understood how to 
obtain the geometry associated to a given Af = 1 gauge theory. In particular, one 
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engineers the gauge theory on D-branes in a Calabi-Yau geometry. Via a geometric 
transition the system has a dual description in terms of a deformed Calabi-Yau 
geometry with fluxes and no D-branes and this geometry is captured by a certain 
Riemann surface encoded in a specific matrix model. 

The matrix model associated to our theory is defined by the partition function 

Z = -—^- [[dX} m [dY}^[dZ}^ exp[--Tr(XYZ -qXZY + V(X))} . (4.1) 
\ol[U(N)\ J 9s 

where there is matrix for each chiral field of the theory. As usual the dimensions of 
the matrices indicated by N are not the same as the fields in the parent field theory. 

As in earlier works [7,20,21] and references therein, the holomorphic integral is to 
be thought of as a contour integral and in practice is performed by choosing Y = 
and then integrating them out. The resulting one- matrix model can be re- written as 
an ordinary A-dimensional integral in terms of the eigenvalues {x^} of X: 

Z= [Udx, } [ ' r r <- r J c _i Ei (42) 

J tl UijXi-qXj 
4.1 Saddle point condition 

The connection with the holomorphic sector of the gauge theory and the associated 
Riemann surface emerges from the so-called genus zero contributions to the matrix 
integral when evaluated in a 't Hooft large- N expansion around a saddle-point. In 
this limit, N — ► oo and g s — > with S = g s N fixed. The saddle point configuration 
of the eigenvalues {x{\ is determined by extremizing the effective action fl4.2| ) leading 
to the condition 

2 y — — y — - — y — = v( Xi ) . (4.3) 

3+i J j J j J 



9s 



In practice, one starts with a classical saddle point solution (of the field theory), 
where the left-hand side of ( |4.3|) vanishes, described by a set of eigenvalues {x a } with 
a certain degeneracy N a , ^2 a N a = N. In the matrix model, one replaces iV a — > N a , 
J2 a N a = N, and takes A\ — > oo and g s — > with each S a = g s N a fixed. In this limit, 
the basic assumption is that the eigenvalues at x a spread out in a one-dimensional 
distribution on an open contour C a in the complex x-plane. This distribution can be 
described by a unit-normalized density p(x) with support on the union of contours 
C = U a C a . The density of eigenvalues can be encoded in the resolvent function 

u(x) = lim iV = [ dy^L . (4.4) 

N->ao N . x- Xi J c x-y 
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The resolvent u(x) is analytic and only has cuts along the equilibrium distribution 
with a discontinuity that determines the density: 



uj{x + e) — uj{x — e) = 2nip(x) , x G C , (4.5) 

where e is a suitable infinitesimal. Integrating the resolvent over any one of its 
branch cuts C a gives us the fraction of the total number of eigenvalues, N a /N , that 
are distributed along that cut. 



The saddle-point equation ( |4.3| ) in the large N limit can now be written in terms 
of the resolvent function as 

S [uj(x + e) + lo{x — e) — quj(qx) — q^ 1 uo{q^ 1 x)] = V'{x) , x G C , (4.6) 

where e is a suitable infinitesimal so that x ± e lie above and below the cut at x G C. 
The content of this saddle point equation is easily understood when recast in terms 
of a new function t(x) (see also [7,20-22]) which naturally encodes the geometry of 
the quantum deformed moduli space and will also allow us to see the duality cascade 
of the field theory in a simple way. The form of the saddle point equation (|4.6j ) 
motivates us to define 

t(x) := f(x) + Sx (q^Luiq^x) - u(x)) , (4.7) 



where f(x) was introduced in (|2.14|) . From the analytic structure of the resolvent 



io{x) it follows that t(x) has cuts at C = U a C a and its rotation by q: q C = U a (q C a ) 



The saddle-point equation (fOl ) is then very simple: 



t(x ± e) = t(q(x =p e)) , xeC. (4.8) 

This is simply a gluing condition that glues a point x ± e on the top/bottom of the 
cut C a to a point q(x =F e) on the bottom/top of the cut q C a . So t defines a Riemann 
surface S which is a copy of the complex x-plane with the cuts identified as above. 



Before we describe the solution of the saddle-point equations and the associated 
Riemann surface S it is useful to develop an intuitive picture anticipating some of 
the key features of the function t(x). We first recall from the results of the previous 
section that in a field theory vacuum preserving a U(N + ) x U(N~) gauge symmetry, 
the field X has the classical eigenvalues e M (l, q, . . . , q p ~ l ) with a degeneracy of iV + 
for each, and e M (g p , q p+1 , . . . , q n ~ l ) with a degeneracy iV~. Correspondingly, in the 
classical limit g s —>■ the matrix model resolvent will have simple poles at x G 
{ e ^5 Q e tn ■ ■ ■ 5 Q p ~ le fi} with residue N + /N, and at x G {q p e )M , g p+1 e M , . . . g n_1 e A1 } with 
residue N~/N. This means that in the classical limit the function t(x) defined in 
( |4.7|) has only two poles, at x = e M and x = e^q p with residues g s (N + — N~)e fl and 
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g s {N + — N~)q p e^ respectively. Thus, in this simplified classical limit, we can see 
that the matrix model depends only on one "modulus" g s {N + — N~). The nontrivial 
result which will be argued below is that when the matrix model interactions are 
turned on (g s ^ 0), each of these two poles will get smeared into a single branch cut 
of t(x) and the saddle point condition glues these two cuts together. Importantly, 
the resulting Riemann surface E only depends on the modulus g s (N + — N~). This 
modulus is naturally interpreted as the gluino condensate in the U(\N + — N~\) gauge 
theory on fractional branes at the end of the duality cascade in the field theory. A 
similar matrix model interpretation was found in [23] for the duality cascade of the 
A\ quiver theory. 



4.2 The general solution 



Let us now turn to the solution of the saddle point equation which will determine the 
Riemann surface E. We will keep the discussion in this section general and focus on 



generic deformations V(x) satisfying the requirement (|2.1| ) for generic n with q n = 1. 
Later, we will make all of this explicit for the cases n = 2 with a cubic deformation, 
and n > 2 with a quadratic mass deformation. 



^From the saddle point equation (|4.8|) it is clear that x is a multi- valued function 
on E since it jumps by a phase factor q = e 27 ™/ n at points where the cuts are identified. 
However, this means that u = x n is single- valued on E. So the saddle-point solution 
of the matrix model has led us to a Riemann surface which admits two meromorphic 
functions t and u. Both t and u have poles at u = oo and zeros at u = 0. Furthermore, 
since E is covered once by the x-plane with cuts, it is an n-fold cover of the complex 
u-plane with an n-fold branch cut joining all the sheets and running between u = 
and u — oo. There are additional cuts C a , the image of each C a , which join adjacent 
pairs of sheets in the cover. 

There is a theorem that if t and u are two meromorphic functions on a Riemann 
surface then there exists a polynomial function in two variables such that 9 

F(t,u)=0. (4.9) 

In the classical limit, where the cuts C a collapse into poles we must recover the 
singular classical curve E ( |2.22|) , F(t,u) = 0. It is useful to write 

F(t, u) = F(t, u) + 5(t, u) , (4.10) 

where 5(t, u) is the quantum correction. This correction is constrained by the be- 
haviour of i at it = and u = oo. Note that due to the branch points at both these 

9 For example, in Farkas and Kra [28] Proposition IV. 11. 6. 
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points a good coordinate is x rather than u itself. Since ui{x) = 0(l/x) for large x, 
we have 

t = f(x)+0(l/x). (4.11) 
Correspondingly for small x we have 

t = O(x). (4.12) 

x— >0 

These two conditions constrain S(t,u) to have the form 

5(t,u)= Yl ^"U a t b (4.13) 

6>0,a>l 
na+tb<(n-l)e 

The parameters {jab} are moduli for the solution of the saddle-point equation. One 
can count the number of these moduli, they number precisely the number of singular- 
ities g of X in ( |2.29| ); this is also the (generic) genus of the curve E. In other words, 



there is a modulus for each singularity of the classical curve E. The picture is now 
clear. A general solution to the saddle-point equation involves a deformation of the 
classical curve where each singularity at which pairs of the n copies of the complex 
u-plane touch are resolved into branch cuts. The extent of each branch cut is then 
controlled by a single modulus. This is precisely the quantum geometry advocated 
in [6]. 



Notice that each singularity is only associated to a single cut in the deformed 
geometry. This means that the matrix model geometry can only accommodate a 
single density for each pair of fractional branes 9l + and 9l~. Otherwise one would 
expect, in general, two cuts in the vicinity of each singularity. This conclusion 
matches precisely the physics of the RG cascade. The basic point is the following. 
Suppose that in the vacuum of the field theory iV + > N~ so there is a net number 
of 91 + fractional branes. In this case, the low energy theory after the cascade should 
involve a U(N+ — N~) gauge group which confines. The matrix model has a single 
density for the net number of 9l + branes and, in particular, there is only one glueball 
field. Conversely, if iV + < N~ then after the RG cascade a U (N~ — iV + ) gauge group 
is left which confines. But once again there is only a single glueball field. 

So the message from the matrix model is in our vacuum we should only take the 
net number of fractional branes 

N, ee N+ - N~ (4.14) 

to infinity in the matrix model. If iV M > (< 0) then these represent 91+ (9l~) branes. 
The remaining min(iV + , N~) fractional branes come as 9t+ © m- ~ JK(*) [u^O, 0] and 
can move off into the bulk of X. As advocated in [5,6], these branes should be treated 
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as probes in the matrix model calculation. This makes sense since these degrees-of- 
freedom have a moduli space and so we cannot integrate out the fluctuations around 
them: they remain as probes which feel the quantum geometry deformed by the \N+\ 
fractional branes. The important conclusion from [5, 6] is is that the dynamics of 
the matrix model can be considered in isolation: there is no "back-reaction" from 
the probes. The effect of the fractional branes is to change the geometry felt by the 
branes in the bulk. 

Given that there will be only one type of fractional brane in the matrix model 
from each pair depending upon whether iV+ ^ N~ , we still have to interpret the 
single cut joining the two sheets in terms of a density of eigenvalues. Recall that in the 
classical limit the singularity in question is located at u = w M . After the deformation, 
the singularity opens up into a cut whose image in the cut x-plane involves two cuts 
C M and q p ^C ll . From this we can reconstruct the density for the fractional. Classically, 
the eigenvalues are g- 7_1 e At , j = 1, . . . ,p^ for 91+ and j — + 1, . . . , n for In 
the matrix model these must spread out onto p M and n — p^ cuts, respectively. Let 
us suppose iV+ > N~ then in order to explain the analytic structure of t it must 
be that the p^ cuts are rotational translates of the form q^~ x C^ j = 1, . . . ,p^. In 
addition, the densities along these cuts as measured by the discontinuities in u(x) 
are all equal. This is natural since these cuts are associated to 9^ whose eigenvalues 
are q-'~ l e fl , j = 1, . . . This means that t defined in ( |4.7[ ) is left with two cuts at 
C M and q VtL C ili which are then glued together by the saddle-point equation ( [4.8|) . The 
whole argument can repeated in terms of fractional branes if iV+ < N~ . 

It follows that if is a contour surrounding the cut C M on the curve £ then 

S^g s N, = g s (N+-N-) = -^-<f *^ . (4.15) 

Note that is positive (negative) for iV+ > iV~ (iV+ < N~). The quantities {S^} 
and {jab} form two bases for the moduli space of the curve. The quantity is 
interpreted as the glueball field of the C7 ( | | ) gauge group factor that confines in 
the IR. 

The final problem remaining is to fix the moduli of the curve in terms of the 
underlying coupling of the theory. This is done by extremizing the glueball superpo- 
tential. The form of this superpotential can be deduced from earlier work [7,21,29]. 
Recall that the glueball superpotential can only depend on one glueball field for each 
fractional brane pair. This motivates the expression 

9 8T 9 

u=i p- a=i 
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where we have defined 

IV N+ >N~ 

Pa = { (4.17) 

Here, JF is the first term, "genus zero", in the expansion of the free-energy of the 
matrix model: \nZ = Y^=o-F g g 2 s 9 ~ 2 ■ We have already seen that 5* M can be expressed 
as an integral around the contour An which surrounds the cut near u = u„. The first 



term in Q4.16Q can also be expressed in terms of a contour integral on E. It is the 
variation of the genus zero free-energy (times gj 1 ) upon transporting \p^\ eigenvalues 
in from infinity to each of the Ip^l cuts of u(x) in the x-plane. This can be expressed 
as the difference of contour integral of tdu/u along a contour that runs out from 
infinity along the upper sheet of the cut minus one that starts on the lower sheet 
of the cut C M out to infinity. These two open contours clearly join to form a closed 
contour which can deformed away from infinity. This contour, which we denote 23 M , 
is the conjugate cycle to A^ on E that goes down the cut C M back along the lower 
sheet, up through the branch cut joining to oo and finally out on the upper sheet to 
join the starting point. In fact {A^, 23^} form a basis for the one-cycles of E. Being 
careful with the normalization, 



dTn if tdu 



(4.18) 



and the glueball superpotential is 



% / f tdu f tdu 



n A — * V J-n U J a u 

fM=l 



Now we consider extremizing W with respect to S^. ^From the definition of 
in ( [4.15| ), the one- form 

1 d /tdu\ 



(taken at constant u) has the property 



w v = <V . (4.21) 



Hence, {w^} is a normalized basis of differentials of the third kind on E. The critical 
equations can then be written in terms of the period matrix t^ v of E: 

9 n 9 

Since E has g moduli these g equations are enough to determine E exactly. 



23 



The equations ( |4.22| ) are precisely the condition that £ is an Yu%=i PuN^-fold 
cover of a torus of complex structure r. In order to prove this we need to find a map 
from £ to the torus E T that covers the latter Yl^iP^^ times. For a point p G £, 
the map is 

/p 9 
^2 mod 2-ni, 2vrir , (4.23) 

where p is an arbitrary base point. In particular, the glueball superpotential at the 
critical point is equal to 

i ( f , / tdu f , / tdu\ 



f(x)z dx 



= — 2ix Res 

x 

In the above, we applied a Riemann bilinear relation and used the fact that dz = 
'Yli fi N li 0J IM and that t can be replaced with f(x) in the vicinity of its pole at u = oo 
(x = oo). 

One remaining issue is to include the bulk branes as probes in the matrix model 
calculation and to work out the effect of the fractional branes on these probes. We 
will simply quote the results of [6] here. The classical moduli space of a bulk brane X 
is modified by simply replacing the function f(x) by t(x). The second term in (|4.7| ) 
represents the effect of the fractional branes on the geometry. In other words we 
replace in (|2.20| ) F(t,u) by the deformation F(t,u) = F(t,u) +5(t,u). This defines 



a quantum deformed moduli space space X: 

uvw = (-l) n+1 F(t,u) . (4.25) 



4.3 The case q = —1 with cubic deformation 

As we have noted earlier the theory with q — — 1 and cubic superpotential (|2.3j ) 
exhibits a duality cascade in a certain regime of parameter space. This example will 
allow us to explicitly see in a simple setting, many of the general features argued 
above. There is one pair of fractional branes so that the deformed curve S is a torus. 
Its form will be determined by a deformation of the classical curve (following from 

F(t,u) = t 2 - -u{u-a 2 ) 2 = . (4.26) 

In this particular case, we can also derive the quantum deformed curve directly 
from the matrix model loop equations and since this is interesting we pause to do this. 
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The loop equations can be derived by noting that the matrix integral is invariant 
under generic reparameterizations. In particular consider a variable shift Sxi ~ 
l/(x — Xi) which should leave invariant the integral Eq. ( |4.2|) . Following the general 
algebraic manipulations shown in Appendix A we are led to the following condition 
(for any finite N): 

(t(x) 2 ) = (\x 2 V'{xf + g s Tx 2 (V'( Xi ) - V'(x)) \— i-1 ) . (4.27) 

\ LX — Xi X + Xi-i/ 

i 

In the large- A limit, when all averages are dominated by the saddle-point config- 
uration of the matrix integral, we can treat this equation as an ordinary algebraic 
equation for t. The first point to note here is that, when V(x) is an odd polynomial, 
t(x) is completely determined by a polynomial function on the right hand side of 
Eq. ( |4.27| ). (If V(x) includes even powers, this is no longer true; the right hand 
side is non-analytic and one also needs to consider higher loop equations.) But we 
are interested precisely in the cubic polynomial deformations Eq. ( f2.3| ), and in the 
large- A limit 

t 2 = - e 2 x 2 [x 2 - a 2 ) 2 + 2 x 2 e / p{x) dx . (4.28) 



c 



4 

where on the right hand side we encounter a constant (times x 2 ) which is the first 
moment of the eigenvalue distribution. Now we can clearly see how this is related to 
the deformed (classical) geometry Eq. ( [2.12|) which was the classical moduli space of 



the field theory. Identifying the variable u in Eq. ( [2.12|) with x 2 in Eq. (|4.28| ), the 
above equation reduces to 

F(t,u) = t 2 - ju{u - a 2 ) 2 + -fu = . (4.29) 

where 7 is the constant determined by the first moment of the eigenvalue distribution. 
This is precisely the quantum deformed curve ( |4.10| ). The effect of the deformation 
is to split the double zero at u = a 2 . Thus the function t lives on the two-sheeted 
complex w-plane with two branch cuts. The associated Riemann surface is a torus. It 
now remains to show how the deformation 7 will depend on gauge theory parameters. 
This can be obtained directly by expressing the hyperelliptic curve ( |4.29|) above in 
Weierstrass form as we show below. 



It is also instructive to see how all this follows from the saddle point equation 
( f4.8| ) for t(x). In this case t(x) has two branch cuts, centred at x = ±a and is defined 

as 

t(x) = —-e x(x 2 — a 2 ) — S[xlu(x) + xu(— x)\. (4.30) 

This in turn implies that t(x) = —t(—x) and hence the two branch cuts are images 
of each other. The saddle point equation provides a gluing condition on these two 
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X= 2L2. 2®, 

00 

1 

Figure 1: The branch cuts of t{x) and the gluing conditions on them lead to a parame- 
terization in terms of a torus Ef with complex structure parameter f . 

cuts as illustrated in Fig. 1. The gluing conditions on the two cuts imply that we 
should think of the Riemann surface as the complex plane with a handle which is a 
torus. To see this torus explicitly and connect it to gauge theory parameters we first 
look for an elliptic or Weierstrass parameterization of the problem. (This is exactly 
the problem that was solved in [7].) In this parameterization the torus of complex 
structure f, which we denote Ef, is realized as the fundamental parallelogram in 
the complex z-plane with sides 2u\ and 2uj 2 with opposite identification. Its complex 
structure is f = u 2 /u 1 . There is a natural map x(z) from the torus E f to the two-cut 
complex x-plane, which must satisfy the conditions: 

A — cycle shift : z — > z + 2ui : x(z) — > x(z); t(x(z)) — > t(x(z)) , (4.31a) 
B — cycle shift : z — > z + 2ui 2 : — > — t(x(^)) — > t(x(^)) . (4.31b) 

Hence x(z) is quasi-elliptic, while t(x(z)) is an elliptic function on the torus Ef. Both 
these functions can be determined precisely in terms of the torus variables, from 
the (quasi)periodicity conditions above and the large-x asymptotics determined by 
Eq. ()Oq) . We find, 

x(z) = V p(z) - ei (f) ; * = -| A 3 ; A = [ei(f)]-3 . (4.32) 

Here we have used standard notation, ej(f) = p(u>i) is the Weierstrass function 
evaluated at one of the three half-periods uji,u>2 and u 3 = uj\ + u> 2 - The constants 
of proportionality are uniquely determined by the large-x asymptotics required by 
Eq. ( |4.30| ), which maps to the behaviour in the vicinity of the pole at z — 0. Now, 
since 

u = x 2 (z) = A 2 (p(z) - e 1 ) , (4.33) 

as anticipated in our general arguments in the previous section, u and t will satisfy 
an algebraic relation which in this case follows from the basic differential equation 
satisfied by the Weierstrass function, 

p'\z) = 4(p(z) - e 1 )(p(z) - e 2 )(p(z) - e 3 ) . (4.34) 
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This leads precisely to the algebraic equation ( (4.7|) with 



7 
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e 3 {r) -e 2 {T) 
ei(f) 



(4.35) 



Indeed, we could have deduced this directly from Eq. (|4.29|) by going to an elliptic 
parameterization, but it is useful to see how the loop equation approach matches up 
with the saddle point approach. 

The relation between the modular parameter f of the torus to gauge theory 
parameters emerges upon extremizing the glueball superpotential of the theory. The 
glueball superpotential depends on the following period integrals ( [4.15|) which can 
be explicitly evaluated in terms of the modular parameter f 

s = _J_It± = ±^P (4 . 36) 



and 



OS 

where 




h(f) (4.37) 



(4.38) 



As we saw from our general arguments, extremizing the glueball superpotential 
Eq. (|4.16|) , written explicitly in terms of the expressions above, fixes the complex 
structure parameter of the torus in terms of the gauge theory coupling in a simple 

way 



This is what we obtain from ( |4.22| ) by noting that the single element of the period 
matrix of a torus is its complex structure parameter. 

Finally we have the critical value of the glueball superpotential, which can be 
determined explicitly using the above expressions or directly from Eq. ( |4.24j ) 




^=1^-^-1^,/—. (4.40) 

where f = t/\N + — N~\. The superpotential has a nontrivial expansion in (fractional) 
instantons when Im(f) ^> 1 and vanishes when N + = N~ . This is in line with our 
expectation that in the absence of fractional branes the low-energy theory is the 
Klebanov-Witten conifold SCFT. Note also that the deformation parameter 7 also 
vanishes in this case (this follows from the fact that both e2(f) and 63(f) approach 
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Figure 2: The one-cycles on the u-plane which lift to the compact 3-cycles of the Calabi- 
Yau geometry. 



the same constant value in the limit Im(f) — > oo). An important feature of the IR 
gauge theory is reflected in this critical superpotential, namely that it is not pure 
M = 1 SYM. It has massive modes which are not decoupled from the Af = 1 gauge 
multiplet. The theory does reduce to the pure Af = 1 gauge theory in a decoupling 
limit t — > ioo keeping fixed the dynamical scale ae 1//3 exp(27rzf/3). This is also 
apparent in the structure of the deformed geometry. 

The full deformed Calabi-Yau geometry (Fig. 2) is obtained from the matrix 
model Riemann surface adding uvw to the right hand side of Eq. ( |4.29 ). In this 
deformed Calabi-Yau geometry: 

t 2 = -uvw + —u(u - a 2 ) 2 + wy (4.41) 

the isolated conifold singularity at u = a 2 has been replaced by a 3-sphere of a finite 
size controlled by 7 which in turn is determined in terms of gauge theory parameters 
by Eqs. (|4.35|) and ( |4.39| ). This is the deformation seen in the geometric dual at the 
endpoint of the cascade. 

It is also interesting to compare and contrast the loop equation (the matrix model 
Riemann surface) for the A\ quiver theory and that for our theory with q = —1. 
Recalling the result of [23], the infrared physics at the endpoint of the cascade in the 
Ai quiver theory is controlled by a matrix model curve: 

y 2 = x 2 + fi 2 , (4.42) 

which is the complex x-plane with a single branch cut. There are two one-cycles, one 
of which is compact (the contour enclosing the cut) while its dual cycle is noncom- 
pact (running from a branch point to infinity). This lifts to the deformed conifold 
geometry y 2 = uv + x 2 + fi with one three-cycle of finite size. This is, of course, what 
one expects for the pure Af = 1 SUSY Yang-Mills theory in the IR. 
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In contrast, as we have seen in Eq. (|4.29|) the theory under investigation leads to 
the two-cut complex w-plane with two compact cycles. One of these corresponds to 
the contour enclosing the cut between u = a 2 + 2i^/t and u = a 2 — 2%-J^/e. The 
dual cycle is given by the contour running from u = to u = a 2 + 2i y / 7 y/e, moving 
down to the second sheet and running back to u = (which is the tip of the second 
cut extending from u = to u = oo). In a certain limit (a — > oo) where one of these 
cycles becomes noncompact we recover the physics of the pure M = 1 SYM in the 
infrared. This is also the decoupling limit for the adjoint fields in Eq. (|3.3|). 



4.4 The case n > 2 with quadratic deformation 



This case with quadratic superpotential (|3J|) allows us to embed the duality cascade 
within a relevant deformation of the UV-finite four-dimensional SCFT. We will now 
determine the quantum deformed geometry associated to this theory. Recall that 
classically, the Higgs branch moduli space of this theory is described by the algebraic 
curve Eq. ( |2.20| ), uvw = (—l) n+1 F(t,u). It has [n/2] isolated singularities of the 



conifold type where v = w = and where the (classical) Riemann surface described 
by F(t,u) = becomes singular. Here the classical curve 

11 2i 2 i 

F(t, «) = n (* + - -r^) = ■ ( 4 - 43 ) 

Note that this curve is only a function of u = x 11 since it is invariant under x — > xq. 
The singularities occur at the points u = u^ p \ p < n given by Eq. (|3.7|) where the 
classical curve above has a double zero. Fractional branes can be located at each of 
these conifold singularities. In all there are [n/2] distinct types of fractional brane 
pairs. The effect of sticking fractional branes at a conifold singularity is to split the 
double zeroes, resulting in a branch cut that joins a pair of n sheets of the Riemann 
surface associated to t as a function of u. This is exactly what we encountered in 
the previous example with n = 2. This deformation leads us to the quantum curve 
(1.10) which takes the form 



,2i^2 



[(n-3)/2] 



F(t,u) = H(t + ef^- 2 -ea^-)+u £ 7 X = . (4.44) 

8=1 ^ ^ (1=0 

The 7 a are the quantum deformation parameters which we would like to determine 
using the matrix model. They are implicitly functions of the glueball condensate of 
the field theory. The curve is actually hyper-elliptic since it has the general form 

[«/2] 

F(t, u)=t n + (3^) u 2 + u s a t a = , (4.45) 



a=0 
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where the s a depend on the moduli. In fact, s a is equal to j a , for a = 0, ... , [(n— 3)/2], 
up to a constant shift, while s a , for a = [(n — l)/2], . . . , [n/2] only depend on the 
parameters in the potential. 

Generically when all types of fractional branes are present the curve has genus 
[n/2], as each of the singularities gets opened up into a branch cut. However, let us 
consider the simplest vacuum, which was described in Section 3.2 where fractional 
branes are present at only one of the conifold singularities. This means that only 
one cut opens up leaving a genus one Riemann surface, or torus. In this case, we can 
solve for t and u (and x) in terms of the Weierstrass parameterization. 



The parameterization we seek can be deduced from the saddle-point condition 
(O), which says that the two cuts of t(x) which are images under a rotation by q 
are now glued together. This adds a handle to the complex x-plane giving us a torus 
(Fig. 2) with a complex structure modulus f. As explained in Section 4.3, we think 
of this torus E f as the complex 2-plane modded out by lattice translations. Going 
around the A-cycle of this torus transforms i by a phase x(z) —>■ qx(z) while the 
B-cycle shift leaves x(z) invariant. Clearly this means that the u(z) = x n (z) is an 
elliptic function of z (this is true only because q is an n-th root of unity). These 
operations also leave t(x(z)) invariant. Along with the large- a; asymptotics (which 
maps to small- z), this determines the two functions uniquely: 

u{z) = An ' ^e^/^) ,n) ' t{x{z)) = B ' [p{z) ~ ^ 2uj ^ • ( 446 ) 

where the constants A and B are 

A = a M±A. B = ^A.l±i.^\ (4 . 47) 

2 9[(pn/n) ' vr 2 1 - q 9f( P 7i/n) V ' 

The modular parameter associated to all the (quasi)elliptic functions appearing in 
the above expressions is f, the complex structure parameter of the torus Ef. 



u(z) is an elliptic function with an n-th order pole at z = and an n-th order 
zero at z — 2pui/n. By determining the coefficients of the singular (and constant) 
pieces in the Laurent expansion of u around z = 0, we can rewrite u(z) in terms 
of the Weierstrass function and its derivatives using well-known results for elliptic 
functions (Eq. (B.19| )). This straightforward (but lengthy) procedure allows us to 
obtain the algebraic relation between u and t. For example for n = 3 and p = 1 we 
find 



F(t,u) 



t 3 + 



L 1 — g 2 J 



a 3 

+ 



u 



3a 



L(l-g 
(tt/3) 



tu 



(4.48) 



Ll-gJ 0f(7r/3) 



p'(2cJi/3) u = 
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Figure 3: The Riemann surface for t(u) in the case n = 3. The two one-cycles on the 
-u-plane lift to compact 3-cycles of the Calabi-Yau space. 

This quantum curve has all the properties that we anticipated in and below Eq. ( |4.45| ). 
(For the n = 3 case, where [n/2] = 1, there is actually only one conifold singular- 
ity.) In particular there is only one quantum deformation parameter which alters 
the coefficient of the term linear in u. As we see below the classical limit — absence 
of fractional branes — corresponds to taking Im(f) to oo. In this limit the quantum 
deformation that we have computed above indeed approaches the result expected 
from the classical curve Eq. ( |2.20| ). The classical Riemann surface associated to t(u) 
consists of three sheets joined via the branch cut running from u = to u = oo. 
The effect of the quantum deformation is to open up the singularity at u = a 3 into 
a square root branch cut that joins two of the three sheets (Fig. 3). The algebraic 
equation for the deformed Calabi-Yau space is, as before 

uvw = (-l) n+1 F(t,u). (4.49) 



Once again we can relate the complex structure r to the gauge theory coupling 



by extremizing the glueball superpotential as in ( |4.22| ). For a gauge theory vacuum 
where we have A^ + and A^~ fractional branes of the two different types, we find the 
relation 10 

pr \pr/(N+ - N~) N+ > N~ . 



\{n-p)T/(N- - N+) N + <N 
Hence, the limit where the deformation vanishes i.e. N + = N~ is when Im(f ) — ► oo. 



10 This can be explicitly seen by evaluating S and dF^/dS which turn out to be dh(f)/df and 
fdh(f)/df — h(f) where h(f) = 2ea 2 [(l + <?)/(! — q)]9i(pir / n) / 8[(pir / n) . Extremizing the superpo- 



tential ( 4.1 6| ) then leads to the desired relation between f and the gauge theory parameter. 
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Finally from ( |4.24j ) we have the critical value of the glueball superpotential: 



2 l + q Bxipnr/n) 



W * = _ 2 |j\T+ - N-\ea*—^. " 1 (4.51) 

1 — q u^pn/n) 

where the modular parameter of the Jacobi theta functions, f is determined as above. 
This value of the effective superpotential coincides with that of the mass deformation 
of the /3-deformed theory obtained in [7]. This indicates that at least at the level of 
the holomorphic sector, the low-energy theory at the end of the RG cascade is the 
mass deformation of the U(\N + — N~\) /3-deformed theory (with (3 = 2n/n). Pure 
M = 1 SUSY gauge theory only emerges in a decoupling limit for the extra massive 
adjoint fields. 



4.5 Relation to periods of the Calabi-Yau geometry 



It is interesting to explicitly see the relation of the matrix model moduli to periods of 
the Calabi-Yau geometry. We know that a vacuum of the A/" = 1 theory with U(N + ) x 
U(N~) gauge symmetry is obtained by taking iV + fractional D3-branes of one type 
and N~ of the second type and placing them at a conifold singularity in the (classical) 
geometry specified by Eq. (|2.20|) . The general picture is that gaugino condensation in 
the infrared theory leads to a geometric dual where the S 2 is replaced by a deformed 
S* 3 without any branes. For the B-model topological string in this deformed geometry, 
the genus zero prepotential term is determined by period integrals of the holomorphic 
(3, 0) form. These period integrals of the deformed Calabi-Yau geometry should 
precisely be the matrix model integrals of the meromorphic one-form over the A and 
B cycles of the Riemann surface. We show that this is indeed the case. 

First, in the simplest case with q = —1, consider the holomorphic (3,0) form 



associated to the geometry (|2.12|) , 



e 2 

ft = dt' A du A dv A dw 5( t' 2 + uvw - —u{u - a 2 ) 2 + ju) 

du A dv A dw (4.52) 



2-y u [—vw + ^-(w — a 2 ) 2 + 7] 

The geometry has compact 3-cycles over which the period integrals can be evaluated 
as follows. First one picks an appropriate reality condition on the coordinates (e.g. 
v = w* assuming that a and 7 are real, positive). The 3-cycles can be identified 
by treating v and w as coordinates on an S* 2 which is fibred over an appropriately 
chosen interval on the w-plane. In this way the period integrals over the compact 
3-cycles in the Calabi-Yau, reduce to contour integrals over the branch cuts (Fig. 3) 
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of a meromorphic 1-form in the w-plane, 



2 



du _ f du 



Q = 2vr (p \ -(u-a 2 ) 2 + 1 ^ = 2>K(b t—, (4.53) 
which are precisely the moduli of the matrix model solution discussed above. 



More generally, we can present an argument along the lines of [30]. For the 
non-compact Calabi-Yau given as the hypersurface in C 4 , 



uvw = {-l) n+1 F(t,u) (4.54) 



where F(t,u) is defined as in Eq.( 4.10 ), we think of the Calabi-Yau space X as a 
fibration over the (t, w)-plane, where the fibre is the curve uvw = (— l) n+1 F. The 
holomorphic (3, 0) form can be chosen to be 

^ dt A du A dv , , 

Q = . (4.55) 

uv 

This has periods over the three-cycles in X which reduce (by Cauchy's theorem) to 

dt Adu . . 

(4.56) 



D U 

where D is a real two-dimensional domain in the complex (t, w)-plane such that 
dD C S. Here S is the Riemann surface F(t,u) = on which the fibre degenerates. 
By Stokes' theorem these integrals reduce to 



f (4.57) 



integrals of a meromorphic one-form over the one-cycles of the Riemann surface S. 
These are precisely the moduli (S^dTo/dSfj) of the matrix model (Eqs. (|4~T5|) and 

CH-) 



5. A Generalization 



Finally, we turn to a more general class of renormalizable deformations of the super- 
conformal field theory (with q n = 1) described by the superpotential 

W = Tr[XYZ - qXZY - m x YZ - m 2 XZ - m 3 XY] . (5.1) 

By shifting each field by a multiple of the identity one has the equivalent description 

W = Tr[XYZ - qXZY - 2aX - 2(Y - 2r]Z] , (5.2) 



33 



where 

m 2 m 3 rxixm-z m x m 2 lr 

This is precisely the deformation considered in [15,16]. It completely resolves the 
orbifold fixed lines leaving behind (n — 1) isolated conifold singularities. 

As previously it is useful to consider a slightly more general form of the theory 
by introducing a general superpotential for X: 

W = Ti[XYZ - qXZY - 2(Y - 2r]Z - V{X)} , (5.4) 

We are going to show that this more general deformation can easily be incorporated 
into our existing formalism by a simple trick. The idea is to "complete the square" , 
in Y and Z by defining 

Y = Y — ^-X- x , Z = Z- -^X- 1 . (5.5) 
q — 1 q — 1 

This move changes the superpotential for X to 



t 



v{x) = ~ x3 - -^j^ 1 ■ ( 5 - 6 ) 

mod n ^ ^ 



Our previous formalism now applies with V replaced by V. Hence the classical 
moduli space X of the Higgs branch is now 

uvw = (-l) n+1 Fy(t,u) , (5.7) 



where Fy(t,u) is given as in (|2.20| ) but with f(x) defined as in fl2.14j ) with V(x) 
replaced by V(x) and 

v = -TiY n , w = -TrZ n . (5.8) 
n n 

Since 

„ = «,+ ( *L)"I, „ = fi+ («L)-I, (5.9) 
VI — qJ u VI — qJ u 



we can write the moduli space in terms of the original variables as 

«» = (-mt, u) + + ( T ^)"u, - . (5.10) 



It turns out that the right-hand side is a polynomial in t and w since the 0(l/u) 
terms cancel. 

The classical singularity structure of X is somewhat different from the case with 
( = rj = 0. The singular line at u = is now resolved and only isolated conifold- 
like singularities remain. In this case the Higgs branch becomes disconnected from 
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the Coulomb branch. As before, the isolated singularities are associated to pairs 
of representations The matrix model analysis can be followed as in Section 4 
to find the quantum deformed geometry. The only difference is that the quantum 
deformation of the curve is slightly more general with ( |4 . 1 0| ) replaced by 

6(t,u)= ^a b u a t b ( 5 ' n ) 

6>0,a>0 
na+eb<(n-l)e 



For example, in the case with q — — 1 and with V(x) = 2ax, the shifted potential 

is 

V(x) = 2ax+ — . (5.12) 

x 

The classical moduli space ( [2.20P is 

uvw = —t 2 + a 2 u + ( 2 v + r/ 2 w — 2a(r/ , (5.13) 

in particular, there is a single singularity at v = w = and u = (t]/<J- The quantum 
curve S is 

F(t, u)=t 2 - a 2 u - CW- + + 7 = 0, (5.14) 

u 

where 7 is the quantum deformation. This curve is, as expected, a torus. It can be 
written in Weierstrass form (|4.34|) by taking 

u = Bp(z) , t = ^4 ' ( 5 - 15 ) 

where 

a = (77^2) 1/2 ' M^f. < 5 - 16 ' 

\( v y v a 93 J 

Here, (72 and are the Weierstrass invariants defined in (p.6j ). The quantum defor- 
mation in ( |5.14j) is then 

1 = 92 (a^> ~ 2aCv - (5 ' 17) 

The critical point of the glueball superpotential will set the complex structure 
of the torus f as in ( Q9|) . Replacing F by F in ( p.!3|) then gives the quantum 



deformed geometry X as a function of the underlying couplings. 



6. Conclusions 



We have studied relevant deformations of the M = 1 SCFT related to M = 4 theory 
by the exactly marginal /3-deformation. We have seen how at intermediate energy 
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scales, certain vacua of the theory (at weak gauge coupling) can have a description 
in terms of an Af = 1 theory undergoing a Klebanov-Strassler duality cascade which 
terminates in an IR theory with a mass gap and gaugino condensation. We under- 
stood all of this rather precisely in the language of the Dijkgraaf-Vafa matrix model 
which allowed us to compute the quantum deformed moduli space seen by a probe 
D-brane at the end of the cascade. Specifically we were able to compute exactly 
the deformation parameters as functions of the microscopic parameters of the gauge 
theory. A natural question arising in the context of these theories, which we hope 
to address in future work, is what is the supergravity/string theory dual for these 
theories. In the UV the field theory (at large- N) is (super)conformal and is dual to 
the near horizon geometry AdS$ x S 5 /T (r being the orbifold action) of D-branes 
at the C 3 /Z n x Z ra orbifold with discrete torsion. 11 The discrete torsion should be 
encoded in the boundary conditions for massless twisted states of the orbifold. As 
noted in [6] the relevant deformations we have studied correspond to turning on cer- 
tain fields in the twisted sector. This deformation should cause the AdS geometry to 
smoothly match on to a Klebanov-Strassler type geometry that describes the cascade 
phenomenon in the IR field theory. It would be extremely interesting to obtain such 
a SUGRA flow. Of course, from the point of view of the field theory results we should 
also expect to find a flow interpolating between AdS 5 x S 5 /T and AdS 5 x X^ 1 ' 1 ) which 
is the conifold theory of [24]. 

Acknowledgements: We would like to thank N. Dorey, G. Policastro and F. 
Quevedo for discussions. 



Appendix A: Matrix model loop equations 



We now outline how the quadratic loop equations follow from the integral ( f4.2| ). 
Performing a variable change 5xi = the invariance of the integral under general 
variable changes leads to the condition 



It- — t. ■ 



n v n^) 

{ S ^ x-xj ' 



x i) 



-y I 



(A.l) 



) 



11 When n is large and q is close to 1, we can also think of this theory as a small deformation of the 
M = 4 theory and consequently its SUGRA dual corresponds to switching on a non-normalizable 
supergravity mode in the AdS§ x S 5 background. The relation between the two descriptions has 
been discussed in [17]. 
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Now, the identities 

1 _ 1 1 

(x — Xi)(xi — q Xj) (x — Xi)(x — q Xj) (x — q Xj){xi — q Xj) 

and 

1 1 q 
= 1 - 



(A.2) 



(A.3) 



can be used recursively along with Eq.( |A.l| ) to obtain the following condition on the 
resolvent function 



n-l „ _ , N n—1 



fco J SN ttt"o X( l- Xi 

We can rewrite this as 

n—1 _ „ n—1 



r 1 n2 1 

[g^(*g fc ) - g fc+1 x^(xg fc+1 ) + ^/(g fc+1 *)J ) = — £ / V +1 *)+ 

fc=0 ' fc=0 

Jy i=l fc=0 ^ 



(A.5) 



In terms of the function t(x) defined in Eq.( |4.7| ) the quadratic loop equation simplifies 



N n-l 



B«v*» = E /v*) + 4 E E ^'( ^"■^ - >• ( A - 6 > 

fc=i fc=i i=i fc=o y ' 1 



For n = 2 and a perturbation that satisfies the condition (|2.3| ) the right hand side 
is a polynomial and the quadratic loop equation is sufficient to determine t(x) in 
the large-iV limit when the loop equation can be treated as an ordinary algebraic 
equation. For n > 2 the loop equation above does not determine t(x). 

For n > 2 we need the higher loop equations. These are more difficult to derive 
directly. However, we can deduce their form by noting an interesting connection to 
the loop equations for the matrix model associated to the A n -\ quiver theory. (Here 
we mean the H = 2 quiver theory with polynomial superpotential deformations for 
the n adjoint chiral multiplets. These softly break M = 2 supersymmetry to M = 1.) 
The higher order loop equations for the A n _\ quiver matrix model were derived 
in [33] using free boson techniques. The quadratic v4 n _i loop equations derived by 
them (see Eq.(2.64) of [33]) is in fact identical to our Eq.( |A.5|) provided we identify 
the k-th resolvent of the quiver model u>k(x) with q k u(q k x) in our theory. We need to 
make a similar identification of the deforming polynomials. Basically, the different 
factors in the quiver are constrained to be the same, up to a phase factor q between 
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adjacent factors in the quiver. This is not a coincidence. It is related to the fact 
that the theory on D-branes at the orbifold C 3 /Z n x Z ra with discrete torsion can 
be obtained by the "orbifolding of an orbifold", [15] (i.e. Z n orbifold of C x C/Z n ). 
The phase factor q between adjacent factors of the quiver represents the single unit 
of discrete torsion. The n-th order loop equations for our theory can be obtained by 
applying this prescription to the results of [33] . It would also be interesting to apply 
these ideas to the derivation of loop equations from anomaly chains as done in [34]. 



Appendix B: Some Properties of Elliptic Functions 



In this appendix we provide some useful — but far from complete — details of ellip- 
tic functions and their near cousins. For definitions and a more complete treatment, 
we refer the reader to one of the textbooks, for example [31]. An elliptic function 
f(u) is a function on the complex plane, periodic in two periods 2ui and 2u;2 • We 
will define the lattice Y = 2uj\L © Iuj-iL and define the basic period parallelogram as 

V = {u = 2/ilu 1 + 2vlu 2 , < n < 1, 0<z/<l}. (B.l) 

The complex structure of the torus defined by identifying the edges of V is 

r = uz/ui (B.2) 

and we also define 

q = e i7TT . (B.3) 



B.l The Weierstrass function 



The archetypal elliptic function is the Weierstrass p(u) function. It is an even func- 
tion which is analytic throughout V, except at u — where it has a double pole: 



p(u) = — + V c k+1 u 2k , 



U4) 



k— 2 

92 93 a j ^ a 

C2 = , C3 = , Cu = — — — > CjCh-j K > 4 

- 20' 28' (2k + !)(£;- 3) ^ 3 3 



3=2 

The Weierstrass function satisfies the fundamental identity 

/ d \ 2 

KduF^J = A ^ u f - 92p(u) - 93 (B.5) 

which defines the Weierstrass invariants #2,3 = 92,3 (wg) associated to the torus. If the 
values of p(u)i) be e i: (i = 1, 2, 3), then {ej} are roots of the equation 4:t 3 —g 2 t—g3 = 0. 
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Here 003 = uj\ + uj 2 . From the formulae connecting roots of equations with their 
coefficents, it follows that 

ei + e 2 + e 3 = 0, e x e 2 + e 2 e 3 + e 3 ei = —92, eie 2 e 3 = ^ 3 . (B.6) 
B.2 The Weierstrass zeta function 



We are also interested in other functions which are only quasi-elliptic. First we have 
((u). It is an odd function with the quasi-elliptic property: 

C(u + 2u e ) =C(«) + 2CM • (B.7) 

Its derivative gives minus the Weierstrass function 

p(u) = -C(u) . (B.8) 

It follows that ((u) has a simple pole at u — 0. Useful identities are 

w 2 C(^i) - ^lCM = y ' ( B - 9a ) 

C(wi + w 2 ) = C(wi) + CM- (B.9b) 

The (quasi)-elliptic functions evaluated at rational multiples of their periods satisfy 
special relations. For example, one useful identity that we need is 



2 



3 

B.3 The Jacobi theta functions 



1/2 



aW3)-^U)= ^(W3) - (B.10) 



We are also interested in the Theta functions 9i(x\r), or 9i(x,q), % = 1,2,3,4. They 
are also quasi-elliptic functions on V when x = nu/2u)i. Each of them satisfies the 
heat equation 

„mp + i mm = ^ ( 

ox or 
They are related to the previous functions; for instance, 



(B.12) 

a;=7TM/2a;i 



o;i 2a>i ^i(x|r) 

where the derivative is with respect to By differentiating this relation one can ob- 
tain further identities relating the Jacobi theta functions to the Weierstrass function 
and its derivatives. 



Another useful identity that we need: 

0;(7t/3|t)-| 2 / tt 



0i(tt/3|t)J 



p(2wi/3). 



(B.13) 
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B.4 "g" -expansions 

To take classical limits (Im(r) — > oo) we use the following expansions 

p( u ) = -^±2l + ( ) cosec 2 2( — ) V - - cos . B.14 



n=l 

Using this and p(u) = — C(u) an d Eq.( |B.12l) , the g-expansion for 9[/9i can be easily 
derived. 

For evaluating classical limits of e^, the following expansion is useful 



OO „ OO ' 

P(u)=(— --+ ^ cosec f—— v \- ^ cosec — - tt . (B.15) 



2^i 



3 ^— ' V 2cji / ^— ' u>i 

n=—oo n= — oo 



B.5 The expression of elliptic functions by means of Weierstrass functions: 

Here we quote the result in article 21.5 of [31] using Eq. flB.12|) to obtain a slightly 
different form. Let f(z) be an elliptic function with a fundamental set (i.e. modulo 
translations by a period) of zeros (a%, a 2 , ■ ■ ■ a n ) and poles (fli, /? 2 , • • • ft n ) so that 

n 

J2(<*r ~ Pr) = 0. (B.16) 

r=l 

Then 

r=l 

where A is a constant. Now if 



Y^Ar^Z-Pr)-™ (B.18) 



m=l 



be the principal part of f(z) at its pole /3 r , then 



m r — 1 



r f \m+l 4 rl" 1 ^ 1 

f(z) = A' + J2{- A r , lC (z { } J;'™' dzm -M * - &)} (B.19) 

r=l m=l 

where A' is an appropriately chosen constant. 
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